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Abstract 

The transfer matrix method is used to analyze resonances in Randall- 
Sundrum models. Although it has successfully been used previously by 
us we provide here a comparison between the numerical and analytical 
models. To reach this we first find new exact solution for the scalar, 
gauge, Kalb-Ramond and q— form fields. Them we calculate numerically 
the resonances by the transfer matrix method and compare with the ana- 
lytical result. For completeness, this is done for models with and without 
the dilaton coupling. The results show a perfect agreement between the 
analytical and numerical methods. 
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1 Introduction 



One of the main models treating aspects of physics of extra dimensions is the 
Randall-Sundrum model [IJ[2] . This model gives a possible solution to the hier- 
archy problem and tell us how gravity is trapped in our four dimensional world. 
However, the lack of a more physical application have posed some challenges 
to this model through the past years. One of its problems is the appearance 
of spacetime singularities due to the presence of infinitely thin membranes. In 
such case, more realistic models based on smooth solutions representing the four- 
dimensional membrane embedded in a higher dimensional spacetime is needed. 
In special cases, where the membranes are generated by scalar fields, its is very 
simple to obtain solutions through the superpotential method. In particular, 
several membrane types were considered: those generated by models with more 
than one scalar fields, deformed membranes and so on [3||9]. Beyond that, the 
analysis of how the localization is made for various fields (tensors and spinors) 
have been understood |To|- 17 . Other aspects include discussions about the 
tensions of membranes |18||19| . 

Another interesting problem is related to understand details of the interac- 
tions between the membrane and several particles that are not the usual zero 
modes. In other words, the approach to this problem considers the computa- 
tional aspects of resonances for various models. In all of the models the massive 
spectrum is determined by a Schrodinger like equation with a potential that 
falls to zero at infinity. The spectrum is not discrete and we have a ill defined 
effective action. Despite of this, just like in the case of semiconductors het- 
erostructures, there is the possibility of appearance of resonances. This analysis 
have been done extensively in the literature [5j[l4, 15, 22 33 . In order to ana- 
lyze resonances, we must compute the transmission coefficient (T), that gives 
a clear physical interpretation of what happens to a free wave interacting with 
the membrane. The idea of the existence of a resonant mode is that, for a given 
mass, the transmitted and reflected oscillatory modes are in phase inside the 
membrane. The transmission coefficient has a peak at that mass value, meaning 
that the amplitude of the wave-function has a maximum value at z — and the 
probability to find this KK mode inside the membrane is higher. This method 
has been used previously to analyze resonant modes of gravity, fermion and 
form fields [34|[35| . 

The background considered here consists of a symmetric Zi thick domain 
wall interpolated between two BPS vacua. The scalar field is important because, 
in some cases, its behaviour is very similar to the gravitational field's behaviour 



37 . The gauge vector field is an important ingredient of the standard model 
and, even being a not localizablc field in theories with conformal symmetries, we 
may understand how its resonances appear due to the specific brane chosen. In 
this case, we consider the corrections due to a dilaton field coupled to the gauge 
field and others antisymmetrical tensors fields. These antisymmetric tensor 
fields arise quite naturally in string theory 38 39 and supergravity [40] and 



play an important role in the dualization processes 41 42 . In particular they 



appear in the R— R sector of each of the type II string theories. These tensor 
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fields couple naturally to higher- dimensional extended objects, the D— branes, 
and are important for their stability. From a more mathematical point of view, 



they are related to the linking number of higher dimensional knots 43 . The rank 



of these antisymmetric tensors is defined by the dimension of the manifold 44 



Beyond this, these kind of fields play an important role in the solution of the 



moduli stabilization problem of string theory 45 ■ 47 



For the physics of extra dimensions is important to study higher rank tensor 
fields in membrane backgrounds. In this context, the antisymmetric tensor fields 
have already been considered in models of extra dimension. Generally speaking, 
the q— forms of highest rank do not have physical relevance because when the 
rank increases the number of gauge freedom increases as well. Such fact can 



be used to cancel the dynamics of the field in the brane 48 . The spectrum 



mass of the two and three-form have been studied in Refs. 49 and 50 in the 



context of five dimensions with codimension one. The coupling between two and 



three- forms with the dilaton was also studied but in a different context 51-53 
An analysis of localization and computations of resonances due to q-forms have 
been previously made by the authors in [54] . 

Here, in this piece of work we show results of particular importance to situa- 
tions where the gravity backgrounds are smooth generalizations of the Randall- 
Sundrum model and procedures that are applied even when the Schrodinger 



potential is not known analytically 20 21 . More precisely, wc continue the 



analysis made in 36 in order to understand the behaviour of scalar, gauge 
vector fields and more general g-forms in a model that can be solved analyti- 
cally. Another goal is to make a comparison between the calculation obtained 
through the solution of the analytical method and the solutions obtained by 
the method of transmission coefficients. The importance of such a model with 
analytical solution resides in the fact that in the models with thick membranes 
their thickness is a constant parameter. When the thickness is parametrized, 
there is a way to obtain the thin limit and the thick membrane limit. In this 
sense, the physics related now can be understood by means of a more general 
situation that can mimics several types of membranes. Through this compari- 
son process we can establish the Transfer Matrix method, once an for all, as a 
tool to understand some physical aspects of extra dimensions. 

The paper is organized as follows. In the first section we review the method of 
computation of transmission coefficients in order to make comparisons between 
different formalisms. In the second section we present the general procedure 
to solve the Schrodinger equation coming from models of extra dimensions and 
solve it for the scenario proposed. In the third section we quickly remember 
the gravitational results and make a generalization regarding solutions as plane 
waves far from the membrane. The remaining sections are devoted to studies 
about the several bosonic fields, i.e., the antisymmetric gauge fields that usu- 
ally appear in literature. They are the scalar field, gauge vector and tensor 
(Kalb-Ramond) fields, besides more general g-forms. Finally we present our 
conclusions and perspectives of work. 
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2 A Simple Example 



In this section we give a review of the Transfer Matrix method and present a 
model with known analytical solution. After that we compare results of both the 
analytical and numerical transmission coefficient computations. The analytical 
model discussed here is already known and can be found in [55] . It is important 
to say that all potentials considered in this manuscript are volcano like, see Fig. 

m 

U(z) 



Figure 1: General potential with parity symmetry with linx^ioo U(z) = 0. 

The potential given in Figjl] can be approximated by a series of potential 
barriers. In each region showed in the Figj2] the Schrodinger equation can be 



2,-1 Zi Zj+1 

Figure 2: The multistep regions, 
solved for the interval < z < Zi, where the potential can be approximated 
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by 

I7(a!) = U{zi- X ) = Ui- X , z 4 _x = ( Zi + *_i)/2. (1) 




Figure 3: The scattering of incident wave into reflected and transmitted waves. 

solution choice is very important for the analytical solution considered in the 
next section. As the final goal here is to compare the analytical and the numeri- 
cal results, we will consider plane waves as boundary conditions to our analytical 
solution. Therefore, the solution in each interval is 

= Ai^e ik ^ z + B i - l e* k *- 1 ', fcj_, - y/X-Ui^, (2) 

and the continuity of the ipi-i{ z ) an d at z — Zi gives us 

Ai ^=MJ (3) 



Bi J \ Bi_x 

In the above equation we have that 

M - JL M fc * + h-i)e- l( - k '- k -^ (h - ki^e-^^-^ \ 
1 2ki \ (h - k i _ 1 )e i ^+ k '-^* (ki + k i - 1 )e i V*- k *-J** J 1 ' 

and performing this procedure iteratively we reach to 

An ^=m( (5) 



Bn J \ Bq 

where, 

M = M N M N _ X ■ ■ ■ M 2 M X . (6) 
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The transmission coefficient is therefore given by 



T= 1/|M 22 | 2 . 



(7) 



In order to numerically obtain resonances we choose z max to satisfy U(z max ) ~ 
10~ 4 and let m 2 runs from U m i n = U(z max ) to U max (the maximum potential 
value). We divide 2z max by 10 4 or 10 5 such that we have 10 4 + 1 or 10 5 + 1 
transfer matrices. 

Now we turn our attention to a model with analytical solution. Consider a 



Schrodinger equation with a potential given by 55 

U(z) 



with profile as in Fig. [4] 



Un 



cosh az 



Uiz) 
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/ 0.02 
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(8) 



Figure 4: The profile of the potential in Ref. 55 with Uq = 1/16 and a = 1 



The analytical solution to the wave function for the potential given by Eq. 
is 

v> = a - er ik/2a F(-i- - s , -»* +8+i, -i- +1,^(1- o), (9) 

where 



^ = tanhaz, k = V2mE/h, 
1 ( 1 , /; 8m^ \ 

s = 2 r 1+ v 1_ w 



(10) 



As stressed in [55] this solution has the desired asymptotic plane wave form. 
After some calculations we can find the transmission coefficient: 
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T = 



T 



sinh 2 (7rfc/a) 
sinh 2 (7rfc/a) + cos 2 [\ <Jl - 8mU a /h 2 a 2 

sinh 2 (7rfc/a) 
sinh 2 (7rfc/a) + cosh 2 '[\^/8mU a /h 2 a 2 - 



8mU /h 2 a 2 < 1, 



-, 8mU /h -V > 1. (11) 



Now we can compare the results obtained by the analytical and numerical 
calculations. We use a = h = m = 1 and Uq = 1/16. Fig. [5] shows: (a) 
The graphic for the analytical transmission coefficient Ta, (b) The graphic for 
the numerical transmission coefficient Tn and (c) The graphic with the ratio 
Tn IT a- There is no distinction between both methods, the agreement is almost 
complete. 




0.01 0.02 0.03 0.04 0.05 0.06 
(b) 
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1.0000 
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0.00 0.01 0.02 0.03 0.04 0.05 0.06 



Figure 5: Transmission coefficients plots in the Landau book example: (a) The 
analytical transmission coefficient Ta, (b) The numerical transmission coefficient 
Ta, and (c) The ratio T N /T A - 



3 The General Analytical Solution 

In this section we show the procedure used through all this work. In models 
describing interactions between fields and membranes, the equations of motion 
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coming from the extra-dimension dependence plays an important role, and have 
a Schrodinger like final form. To arrive at that this result we provide a short 



review here (details can be found in 37 ). First of all we consider the line 
element of the background space as 

ds 2 = e 2A ^i hv + e 2B ^dy\ (12) 

where rj^ = diag(— 1, 1, 1, 1). This metric comes from the inclusion of the 
interaction produced by a dilaton field that, in a specific model, gives rise to 
localization of gauge fields. 

The Einstein's equations give the relation B — A/A between the functions 
A(y) and B{y) [37]. These equations are obtained from the usual superpotential 
method in models containing kink like defects. It is important to cite that the 
dilaton field have been already been treated in the Randall-Sundrum scenario 



53 . As we are going to consider models with and without the dilaton (B = 0), 
we introduced a constant b to take account of this and we get B{y) = (l—b)A(y). 
The value of b with the dilaton is 3/4 and 1 without the dilaton. Using this 
metric background, in general, the equations we need to work have the following 
form 

(-^ + P'^iy + Mv) = ™ 2 Q(vMy), (13) 

where P{y) = r yA{y), Q(y) = e~ 2bA ^ and V(y) = for all fields except gravity, 
in which V(y) — 2A"(y) — 2(1 + b)A'(y) 2 . Here, 7 is a constant that depends 



on the field. We can transform ( 13 ) in a Schrodinger like equation through the 
transformations 

% = Kv)> *Kv) = WvMz), (u) 

with 
and 



f(y) = VQbjj, n(y)^exp(P(y)/2)Q(y)- 1 /\ (15) 



U(z) = V(y)/f + (P'(y)n'(y) - Q"(y)) /nf, (16) 

where the prime means the derivative in y. We emphasize here, that the above 
expression is useful when dz/dy = f(y) do not have an analytical solution. 
When this solution is known it is better to express the potential in terms of 
derivatives in the z coordinate: 

il(z) il{z) f(z) 

where f(y) — f(z). When the above steps are performed we get a Schrodinger 
like equation 

{-^ + U}i,(z)^m^(z), (18) 
with potential U(z) given by 
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U(z)=cA (z) + c 2 {A'(z)) 2 , 



(19) 



where c = — (7 + b)/2 for all fields and c = 3/2 for gravity. 

An analytical solution has been obtained previously for the gravity case 



36 but without the dilaton coupling. We use the same method to solve the 



case with the dilaton coupling and for other hclds. The relation between the z 
and y coordinates in all cases is given by ^ = e' a ' = e~ bA ( z >, for \z\ > d/2, where 
d is a constant. With this, we obtain 



A(y) = A(z) = -- b H\z\+P). 



(20) 



The above expression determines the Schrodinger equation potential in the re- 



gion \z\ > 2 as 



U(z) 



(21) 



with a = I + rj. Since we can have resonance only for positive potential, we 
restrict a to be positive, i. e, c/b > or c/b < —1. For the region \z\ < d/2, we 
choose 



A{z) 



1 



In cos(vVb|z|), 



(22) 



such that 



U(z) = -V , \z\ < d/2. (23) 
Continuity of the metrics and its derivative at z = ±d/2 give us the relations 



Q+p) =cos(v%fr b / c , 

(^+/3) - 1 = lVV t<m(y/V § ). 

In order to obtain Vo, j3 and d we will introduce the parameter x = d^/V^/2 as 
in 



30 



For the region \z\ > d/2, the Schrodinger equation can be solved if we 
consider 'ijj(z) = \f\z\ + j3cj>{z). With this transformation we obtain 



m 2 — 



1 _ a 

M\z\+W ~ (kl+/?) 2 



= 0. 



(24) 



Multiplying now by (\z\ + /3) 3 / 2 and dehning u — m{\z\ + (3) we arrive at 



u 2 4>" + u<t>' + [u 2 -v 2 \4> = 0, 



(25) 
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where the prime means a u derivative and v 2 = (5 + f) 2 - This is a Bessel 
equation of order v. Here we are interested in solutions behaving like plane 
waves when z — > ±00. Solutions with this properties are given by 



'hP(u) =J v (u) + iY v (u), 
Hl 2 \u) =J v {u)-iY v {u), 



with asymptotic behavior 



lim^oo H$p (u) = 
lim^oo (u) = 




(27) 



The choice of asymptotic plane waves is related to the fact that we want 
to compare this solution with the transfer matrix method, where a plane wave 
colliding with the brane is used. It is important to stress here that this is 



different from the solution found in 36 , where an exact solution was used but 



no plane wave was considered. Using this consideration, the solution is given by 

1>{z) = (AHpHu) + BHV\u)) , (28) 
with the desired behavior. 

lim i>(z) = J— (Ae lmz + Be~ imz ) . (29) 
In the region \z\ < d/2, we have the solution 

i)i(z) =ae iKz + be~ iKz , (30) 



with K — v m 2 + Vq . The subscript / means that the solution is in the central 
region and all the above constants must be fixed by the boundary conditions. As 
we arc interested in calculate resonances we must consider a plane wave coming 
from —00. This plane wave will collide with the membrane and will generate 
a reflected and a transmitted wave. Therefore for z < —d/2 we must have a 
linear combination of waves moving to the left and to the right. For z > d/2 we 
must have only one wave moving to the right. We define the solution in both 
regions by iPl(z) and iPr(z) respectively. In order to analyze resonances we fix 
the coefficient of the incoming wave equal to one. Therefore we choose 

feW = J^(H^{u) + B 1 Hi 1 \u)), z<0, 
M*) = \[^A 2 Hl 1 \u) 1 z>0. 

V TO 

Defining E{z) = J^H^ \u),F{z) = J^H^ ] (u) for z < and taking the 



continuity of the wave function and its derivative at z = id/2, we obtain 
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1 



' E{-d/2) F(-d/2)\ 

E'{-d/2) F'(-d/2)J KBiJ ~ \ iKe~ iK i -iKe 




(31) 



and 




(32) 



F(-d/2) 0\ (A 2 \ ( e iK i 
-F'i-d/2) \iKe lK % -iKe 

After some algebra we finally obtain 

-UK 

2 ~ n{2KFF' cos Kd + (F a — F 2 K 2 ) sin Kd) ' (33) 
with the transmission coefficient given by 

2 16K 2 
T= '^ 2 '~ = n 2 \2K FF' cos Kd + (F' 2 -F 2 K 2 ) sin Kd\ 2 ' ^ 

That is the final expression we use as starting point to make comparison 
between analytical and numerical methods. It is important to mention the 
number of iterations made using the method applied here. For calculations 
involving the dilaton contribution we have used 10 5 matrices, 2 x 10 5 matrices 
for the gravity case with x = 1.5, and 10 5 matrices with x = 1.0 for the scalar 
and gauge field cases. We have made another calculation of 2 x 10 5 matrices 
for the parameter x — 0.15 again for the vector gauge field resonances (x is 
the parameter describing the thickness of the membrane) . For the tensor fields 
we computed 10 5 matrices for x — 1.5 and x = 1.0. For the case without the 
dilaton contribution, we have computed 3 5 matrices for x — 1.5 and 10 5 for 
x = 1.0 for the scalar field resonances. In the case of vector fields we made 10 5 
iterations with x — 1.0 plus another computation of 10 5 matrices for x = 0.15. 
It is important to mention that we have found, for the case without dilaton, the 
same results for resonances in the gravity and scalar field cases. The results of 
these calculations will be discussed below. 



4 The Gravity Case 

As mentioned before the gravity case without dilaton has been considered pre- 
viously in [36]. We will reconsider it here for completeness. The potential for 
this case is given by 

Ui(z) = \e 2A ^ (2A"(y) + 5(V(y)) 2 ) • (35) 
The transformation to the z coordinate is obtained using ^ = e~ Al ( z ) to get 

U 1 (z) = ^A l "(z) + l(A' 1 (z)) 2 , (36) 
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what is equivalent to choose 6=1 and c = 3/2 in the previous section. From 
this we can obtain the solution taking v = 2. 

The case with the dilaton coupling has the potential 



U 2 (z) = -e 



= 2. P 3A 2 (y)/2 
2 



M'\v) 



-{A' 2 {y))'> 



Where using the transformation ^ = e 4^2(2) we obtain 
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dy 
3 



U 2 {z) = -A 2 "(z) + 



(37) 



(38) 



what is equivalent to choose 6 = 3/4 and c = 3/2. The solution is given by 
taking v = 5/2. Using now the expression for the transmission coefficient we 
obtain the results with dilaton in Fig. [6j In the figure we can see the appearance 
of one peak of resonance for the parameters x = 1.0 and AV3A/ 3 = 20. These 
peaks of resonance are very dependent on the parameter x, when we take x = 1.5 
and A\/3M 3 = 20 the transmission coefficient for gravity with dilaton presents 
seven peaks of resonances shown in Fig. [7] 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1 



Figure 6: The transmission coefficient for gravity with dilaton for x = 1.0 and 
Av / 3M 3 = 20. The solid line is the analytical calculation and the dashed line is 
the numerical one. 



5 The Scalar Field Case 

Now we analyze the scalar field resonances by the analytical model proposed. 
As in the previous section we must obtain the potential in the z coordinate for 
the cases with and without the dilaton. These are given by 
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0.0 1 m 

0.40 0.45 0.50 0.55 0.60 0.65 



Figure 7: The transmission coefficient for gravity with dilaton for x = 1.5 and 
A V 3M 3 = 20. The solid line is the analytical calculation and the dashed line is 
the numerical one. 



n .< : ) ^ < :U 2 ( "(f + \)Mvf + (~ - l^m 2 ) . (39) 



(40) 



where a = — 15/4— A\/3Af 3 . Performing now the transformations ^ = e 3yl i( z )/ 4 
and g| = e~^ 2 ( z ) respectively, we obtain 

ftW = (-(f + §)4"(*) + (f + f)V.(*)) a 



and 



faW = (s^"W + 7(4(=)) 2 )- ( 42 ) 



Therefore the solutions to both cases can be found by using 

3 ,3 AV3M 3 

h = v <* = (z + —2—), 
3 

h = 1, c 2 = -, 
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and we get 



, _ f 5 2AV3M 3 2 



= <2 + 



V2 = 2. 

We now compare these results with those obtained from the numerical method 
of transmission coefficients. We show in Fig. [8] the transmission coefficient for 
analytical (lined) and numerical (dashed) without dilaton with x = 1.0, showing 
one peak of resonance. In Fig. [9] with dilaton field we still have one peak of 
resonance for the parameters x = 1.0 and AV3M 3 = 20. In Fig 
Transmission coefficient for scalar without dilaton for x 
and 
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in Fig. 



we have the 
1.5 andAV3A/ 3 = 20 

11 showing eight peaks of resonances (the first one near m = 0, can 



be made more visible if we just give a zoom in that region), the Transmission 
coefficient for scalar with dilaton for x — 1.5 and \\J3M 3 = 20 showing three 
peaks. 




0.0^ 
0.8 



0.9 



1.0 



1.1 



1.2 



Figure 8: Transmission coefficient for scalar without dilaton for x = 1.0. Lined 
is analytic and dashed is numeric. 



6 The Gauge Field Case 

Now we turn our attention to analyze the gauge field resonances through the 
analytical context. Again we need to obtain the potential in the z coordinate 
for the cases with and without the dilaton. These potentials are given by 

U l[z) = e 3A l(y) /2 (_q + 3 )Ai „ (y) + { a ^ ){A > i{y)) ^ , (43 ) 
U 2 (y) = e^(^A"(y) + l(A'(y)ry (44) 
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where a = — 7/4— XVSM 3 . Performing now the transformations ^ = e 3j4 i( z )/ 4 
and g| = e - " 42 ^ respectively, we obtain 



f?iW = (-(f + §)*"(«) + (f + g) a (^iW) a 



'(i + ^^"W + (i + ^)*WW)"), (45) 



and 



C7 2 (z)=Qa 2 "(z) + 1(4(z)) 2 ). (46) 
Therefore the solutions to both cases can be found by using 



3 ,1 A%/3Af 3 . 

*i = 4 , c 1 = (- + ^— ), 

&2 = 1, C 2 = -, 



what gives for the solutions with and without the dilaton 



,7 2AV3M 3 2 



// 2 



1. 



We show in Fig. [13] and in Fig. 12 the transmission coefficient for analytical 
(lined) and numerical (dashed) with and without dilaton with x — 1.0 showing 
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0.20 0.22 0.24 0.26 0.28 0.30 0.32 



Figure 10: Transmission coefficient for scalar without dilaton for x — 1.5 and 
X^/SM 3 — 20. Lined is analytic and dashed is numeric. 

one peak of resonance, and in Fig. [14] we show the three peaks of the transmis- 
sion coefficient for vector field with dilaton for x = 1.5 and A\/3A/ 3 = 20. 

7 The Kalb-Ramond Case 

The last bosonic field to be analyzed in five dimensions is the Kalb-Ramond 
field. Now we analyze this field resonances in the same way as in the last 
sections. The potential in the z coordinate for the cases with and without the 
dilaton are given by 

UM = e 3 ^ )/2 (-(f + |W(y) + (f - ^)K(2/)) 2 ) , (47) 

Uv) = z 2M{v) + \(A(y)f) , (48) 

where a = 1/4 — Xy/3M 3 . Performing now the transformations ^ = e - 3j4 i( z )/ 4 
and ^| = e _j4 2(z) respectively, we obtain 

ftW = (-(f + g)*"(») + (f + |) a (^iW) a ) 

= ((-i + ^ 5 )i I "W + (-i + ^(ii(*),(49) 

and 

U2(z) = (-lA 2 "(z) + ±(A' 2 (z)) 2 y (50) 
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10.0 10.5 11.0 11.5 12.0 12.5 13.0 



Figure 11: Transmission coefficient for scalar with clilaton for x — 1.5 and 
X^/SM 3 — 20. Lined is analytic and dashed is numeric. 



Therefore the solutions to both cases can be found by using 



h = 



ci = (- 



i xVup 



bo = 1, 



C2 = 



2' 



), 



which gives for the solutions with and without the dilaton 



ul = (- 

^2 = 



1 2AV3M 3 , 
6 + 3 ' 



The case without the dilaton has a = —1/4 and therefore will not be con- 
sidered because it gives a negative potential. 



We show in Fig. 15 the transmission coefficient for analytical (lined) and 
numerical (dashed) with dilaton with x — 1.0 and in Fig. 16 with x = 1.5. In 
both cases we have AV3M 3 = 20. In Fig. [ 
around m = 11.5 and, in Fig. 16 we have got five peaks. 



we see just one peak of resonance 



8 The q— form Case 

The results of the last three cases can be summarized and generalized to a q- 
form in a p-brane, where p = D — 2. In a recent article the present authors have 



studied the resonances of theses fields numerically in a different background 54 
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0.0 1 m 

0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 



Figure 12: Transmission coefficient for vector field without dilaton for x = 1.0. 
Lined is analytic and dashed is numeric. 



There the potential in the z coordinate for the case with the dilaton is given by 

tj l[z) = e 3A l( ,)/2 + 3 K(y) + ( « _ ^Mf) ■ (51) 



In this equation a = (8q — Ap — 3) /4 — A v / 3M^ . Performing the transformations 

dz = e -3A!(*)/4 we obtain 



ftM - (-(f + |)VW + (f + |)'(*(.))'). (52) 
For the case without the dilaton, we have obtained the potential 

U 2 (z) = e 2 ^> ((g - P -){q - ^)(A' 2 (y)) 2 -(1- f KG/)) , (53) 
and using the transformation |* = e - " 42 ^* 1 we obtain 

U2(z) = ((| - ?)^'(*) + (| - 9 ) 2 K(z)) 2 ) . (54) 
Therefore the solutions for both cases can be found using 

3 ,p AV3M 3 . 

61 = 4' Cl = ( 2~ 9+ ^— } ' 

62 = 1, c 2 = (--q), 

which gives the solutions for the cases with and without the dilaton 
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(-2^2 

1 3 ' 



1+P 



The above formulas can be verified for the cases p = 3 with q = 0, 1, 2 and 
agree with the last sections. For the case with the dilaton we show results for 
the parameter a only. The condition for positivity of the potential is a < —3/4 
or a > 3/4. For the case without the dilaton we only have the parameter 
(p/2 — q + 1/2). For this case the condition for the positivity of the potential is 
given byq<p/2ovq>p/2 + l and therefore we understand why, for the two 
form the potential is negative in p = 3. It is important to note that for p > 3 



higher forms can be analyzed. We show in Fig. 17 the potential profile for a 
g-form field without dilaton with p/2 — q = 5. 

For the g-form case the transmission coefficient is plotted in Fig. [18] for 
x = 1.0 and a — —1 and in Fig. [19] for the same x and a = —2 considering the 
dilaton. For the case without the dilaton the transmission coefficient is shown 
in Fig. [20] for p/2 - q = 5 and in Fig. [2l] for p/2 - q = 7 , x = 1.0 in both 



cases. One can see one peak of resonances in all figures except for Fig. 18 that 
has two peaks. For completeness, we show the transmission coefficient in Fig. 
22 and Fig. 23 without dilaton for p/2 — q = 5 and p/2 — q = 7. They show six 



and eight peaks of resonances respectively. 
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Figure 14: Transmission coefficient for vector field with dilaton for x = 1.5 and 
\y/3M 3 — 20. Lined is analytic and dashed is numeric. 



9 Conclusions and Perspectives 

In this work we have analyzed the transfer matrix method in the light of models 
with analytical solution. The background considered consist of a symmetric 
thick domain wall interpolating between two BPS vacua. It has been shown 
previously that this background allows for an analytical solution of the gravity 



field modes 36 . By following the same lines of reasoning, as a first result we 
have found exact solutions to all the bosonic fields, namely, the scalar, the vector 
and the KR field. More than this, our solution to the Gravity case is slightly 
different of ref. (36 . The reason is that in that work they considered the wave 
function as real. Here we need a complex wave function since we consider plane 
waves arriving from the infinity and colliding with the membrane. 

The first case used to test the numerical method was that of Ref. [55j . For 
this well known case the difference in the transmission coefficient values is of 
the order of 10 -17 . For the other fields namely: the scalar field, the vector field 
, the Kb and q— form fields the agreement was impressively good. The results 
showed that the number of resonance peaks are very sensitive to the parameter 
x and the thickness of the membrane. 

For the gravitational field we show Fig. [6] In that picture we can see the 
appearance of one peak of resonance for the parameters x = 1.0 and A%/3Af 3 = 
20. However, in Fig. [7] we have the Transmission coefficient for gravity with 
dilaton for x = 1.5 and A\/3M 3 = 20 presenting seven peaks of resonances. This 
shows the rule of the x parameter in the setup chosen. We show for the scalar 
field, in Fig. [8j the transmission coefficient for analytical (lined) and numerical 
(dashed) without dilaton with x = 1.0, showing one peak of resonance. In Fig. [9] 
with dilaton field we still have one peak of resonance for the parameters x = 1.0 
and AV3M 3 = 20. In Fig. 10 we have the Transmission coefficient for scalar 
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without dilaton for x — 1.5 and A\/3M 3 = 20 and, in Fig. 11 showing eight 
peaks of resonances (the first one near m = 0, can be made more visible if we 
just give a zoom in that region), the Transmission coefficient for scalar with 
dilaton for x = 1.5 and A V 3A/ 3 = 20 showing three peaks. For the case of 
the vector field, we show in Fig. [13] and in Fig. [12] the transmission coefficient 
for analytical (lined) and numerical (dashed) with and without dilaton with 
x = 1.0 showing one peak of resonance, and in Fig. [14] we show the three 
peaks of the transmission coefficient for vector field with dilaton for x = 1.5 and 
AV3M 3 = 20 

In the case of the Kalb-Ramond field we show in Fig. [l5]the transmission 
coefficient for analytical (lined) and numerical (dashed) with dilaton with x — 
1.0 and in Fig. [l6]with x = 1.5. In both cases we have \V 3A/ 3 = 20. In Fig. 
[15] we see just one peak of resonance around m = 11.5 and, in Fig. [16] we have 
got five peaks. For the g-form case we give the graphics for the transmission 



coefficient in Fig. 18 and in Fig. [Kfl with x = 1.0 to a = —1, —2 with dilaton 



and in Fig. 20 and Fig. 21 with p/2 — q = 5, 7 without dilaton. In the first one 
we can see two peaks of resonances, while in the remaining plots we have found 
just one peak. For completeness, we also give the graphics for the transmission 



coefficient in Fig. 22 and Fig. 23 with p/2 — q = 5,7 without dilaton. They 



show, respectively, six and eight peaks of resonances. 

The transfer matrix method has been extensively used by the present authors 
to analyze resonances. But it was important to compare the numerical studies 
with the analytical ones to show that the numerical method is reliable and have 
a perfect agreement with solvable cases. In here we give this comparison for a lot 
of cases which encloses the presentation of the method. There are still more cases 
to be studied like, for example, splitting membranes, a case which would just 
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Figure 16: Transmission coefficient for tensor field with dilaton for x = 1.5 and 
A V 3M 3 = 20. Lined is analytic and dashed is numeric. 



change the Schrodinger potential. The study of fermionic resonances within this 
formalism is also a perspective of the present work. Another important aspect 
is related to corrections obtained, for example, to the Newton's law. Although 
this is not the main goal of the present work, it is important to mention that 
these corrections can be obtained through the transmission coefficient. The 



issue of metastable states has been better studied in Refs. [20 56 . Generally an 
expression is given in the literature which involves ip m (0), however, as stressed 
in 56 , this is obtained after the spectral density expnasion in a complete set of 
modes. However, for our purpose it is better to keep the correction in terms of 
the spectral density, that is given by 

r e -my 

V(r) ~ / p(m)dm. (55) 

J r 

For a sharply peaked resonance we have p = 5(m — mo) and we recover the 



well known formula for discrete masses 56 . When there is a peak in T, by 
current conservation there is also a peak in the spectral density and the above 
integration can be approximated by a sum. From the knowledge of the spectral 
density other important quantities as the width of the resonance can also be 
obtained. 
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Figure 17: The potential profile for q— form field without dilaton with p/2 — q = 
5. (a) x — 1.0 and (b) x = 1.5. 
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